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Abstract. QCD-inspired phenomenological analysis of experimental moments of 
proton and deuteron structure functions F2 have been presented. The obtained 
results on the d/u ratio at large-x, isospin dependence of higher twists and com- 
parison with Lattice QCD calculations were discussed. We remind shortly these 
results: the obtained ratio is consistent with the asymptotic limit d/u — » at 
X — > 1, the total contribution of higher twists is found to be isospin independent 
and the non-singlet moments are in excellent agreement with the Lattice data. We 
present here some details of the analysis triggered by the public discussion. 

Measurements of the nucleon structure function F2 provide the information 
about the longitudinal momentum distribution of partons. These distributions 
being governed by soft strong interactions cannot be described by perturbative 
QCD methods. Only Lattice QCD simulations allow to evaluate these quanti- 
ties. Recent measurements of proton and deuteron structure function moments 
over wide Q^-interval [1,2] and the evaluation of neutron moments [3] allowed 
to improve the knowledge of these non-perturbative distributions. Detailed 
descriptions of these analyses are given in papers mentioned above, whereas 
in the present proceeding we develop further two arguments selected by the 
public discussion. 

Experimentally extracted moments of the proton and deuteron structure 
functions F2 were analyzed to separate leading twist (LT) and higher twist 
(HT) terms. This was performed by fitting the data Q^-dependence with the 
following expression: 
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where LT„ is the LT part of the n-th moment evaluated at NLL accuracy, 
as is the running coupling constant, /i^ is an arbitrary scale (taken to be 10 
(GeV/c)^), is the matrix element of corresponding QCD operator, 7^ is its 
anomalous dimension, /3o = 11 — f -^f with Np being number of active flavors, 
r is the order of the twist and k is the maximum HT order considered. The 
number of HT terms (fc) in the expansion 1 is, of course, arbitrary because we 
don't know at which 1/Q^ power the series converges. Moreover, anomalous 
dimensions of perturbative cocfhcients in front of HT terms are known in a very 
few cases [5,6] . Most of x-space analyses neglect this dependence assuming 7,^ = 
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for r > 2. In the presented analyses the anomalous dimensions were varied as 

free parameters and extracted from the best fit to the data. The results show 
a very strong sensitivity of the fit to the values of HT anomalous dimensions 
at low-Q^. Indeed, it can be seen in the comparison of two twist expansions 
shown in Fig. 1: one using HT anomalous dimensions as free parameters and 
another one assuming them to be zero. The lower limit of the fitted Q^- 
interval was taken to be 1.2 (GeV/c)^ for the full fit. In the fit with fixed 
anomalous dimensions it was increased to 3.6 (GeV/c)^ by the condition of 
having the same number of degrees of freedom. It is evident that 

only the variation of anomalous dimensions permits to describe the data until 
= 1.2 (GeV/c)^ by two HT terms. This observation emphasizes that the 
knowledge of perturbative anomalous dimensions of HT terms is crucial to 
single out individual HT operator matrix elements. 
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Figure 1: Fit of the structure function moment Ms with Eq. 1 using higher twist anomalous 
dimensions 7^ as free parameters (left) and assuming 7,^ = (right): dashed line - the leading 
twist contribution, dotted lines - twist-4 and twist-6 contributions, dot-dashed line - the total 
higher twist contribution, solid line - the total fit. 



In the extraction of neutron moments we assumed the dominance of the 
Impulse Approximation (lA) in the LT part of deuteron moments and treated 
other nuclear effects as, model dependent, corrections to this approximation. 
This allowed for a simple extraction of LT neutron moments from the following 
algebraic relation: 
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where M^, and are LT moments of the proton, neutron and deuteron, 
respectively. Nj^ is the moment of the nuclear momentum distribution 



i.e. the structure function of the deuteron composed of point-hke nucleons (see 
Rcf. [3] for details). The dominance of lA in LT moments, however, imphes 
that processes beyond lA contribute mainly to HT terms. These processes are 
the scattering off correlated nucleons (Final State Interaction (FSI)) and the 
scattering off a nuclear constituent different than the nucleon (Meson Exchange 
Current (MEC)). Only the lowest n = 2 LT moment, sensitive to the low- a; 
dynamics, carries a small contribution from FSI and MECs estimated to be 
< 0.5 %, whereas it is found to be negligible for higher (n > 2) LT moments. 
In fact, the LT part of FSI shown in Fig. 2 contaminates structure functions 
at very low x {x < 0.1) values because the nucleon spectator has a long time 
^'^ < 1 /Mx (here M is the nucleon mass) to interact with nuclear environment 
while awaiting return of the active quark [7]. However, for higher moments 
(n > 2) the mean x value is close to unity and therefore the time left for the 
interaction ^° < 1/M « l/m,r- Here we assume that nuclear interactions are 
carried mostly by pions. 




Figure 2: Example of FSI mechanism in the inclusive electron-deuteron scattering. The 
nuclear interaction between the nucleon spectator and nucleaj: medium is likely carried by a 

colorless pion. 



Hence, the bulk of FSI was expected to give a contribution to the HT term 
of the moment expansion, where the current quark rescatters from the nuclear 
spectator. In order to test this assumption phenomenologically we used calcu- 
lations of FSI in the quasi-elastic peak region based on the model from Ref. [4] . 
To this end, we computed moments of the modeled deuteron structure function 
F2 including and excluding FSI. The ratios between these two calculations for 
n — 2 and n = 8 are shown in Fig. 3. From the figure it follows that, indeed, 
FSI contribution disappears at large generating an additional HT term, 
while LT part is unaffected by this contribution. Moreover, the size of FSI 
contribution in moments does not exceed few % at lowest analyzed Q^. and it 
is much smaller than the total nucleon HT contribution estimated to be about 
w 25% at the same Q'^. We speculate, therefore, that also the HT term has 
not more than 20% contamination from FSI, and Eq. 2 is also applied to this 



term within w 20% accuracy. This accuracy is comparable to the precision of 
the extracted total HT term [1,2]. 
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Figure 3: Ratio of the dcutcron structure function moments calculated using the parame- 
terization of from Ref. [4] including and excluding FSI in the quasi-elastic channel: the 
solid line - n = 2, the dashed line - n = 8. 



Summarizing, the presented analysis of the experimental moments of proton 
and deutcron structure functions F2 showed that: 

• knowledge of perturbative anomalous dimensions of higher twist terms is 
crucial to single out individual higher twist operator matrix elements; 

• for n > 2 FSI mechanism appears in the nuclear structure function mo- 
ments as an additional higher twist term. Its partial estimates indicate 
that the relative contribution of FSI to the total higher twist term does 
not exceed 20%, comparable to the precision of the total higher twist 
extraction. 
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